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ABSTRACT 


fc-The  theory  of  planar  near-field  measurement  is  discussed  for  two 
different  approaches:  The  scattering-matrix  formulation  approach  and  the 
reciprocity  theorem  approach.  The  interaction  of  test  antenna  and  probe 
antenna  is  calculated,  in  which  is  taken  account  of  the  presence  of  the 
probe  used  to  sample  the  field  distributions  since  the  probe  is  not 
assumed  to  be  ideal  (elementary  magnetic  or  electric  dipole) ,  it  is  said 
that  probe  correction  is  applied. 

It  is  shown  that  both  approaches  lead  to  equal  results.  The  derivation 
of  far-field  antenna  characteristics  from  near-field  data  is  then 
discussed  and  near-future  research  items  are  shortly  mentioned. 
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SAMENVATTING 


De  theorie  van  planaire  nabije  veld  metingen  wordt  besproken  voor  twee 
benaderingswijzen:  De  verstrooiings-matrix-formulering  benadering  en  de 
reciprociteits-theorema  benadering.  De  interactie  van  testantenne  en 
probe  is  berekend,  waarbij  rekening  is  gehouden  met  de  aanwezigheid  van 
de  probe  welke  gebruikt  wordt  om  monsters  van  de  veldverdeling  te  meten 
Daar  de  probe  niet  ideaal  is  verondersteld  (elementaire  magnetische  of 
elektrlsche  dipool),  wordt  wel  gezegd  dat  probe-correctie  is  toegepast. 
Aangetoond  wordt  dat  beide  benaderingswijzen  tot  dezelfde  resultaten 
leiden.  Uiteindelijk  wordt  de  afleiding  van  verre  veld  antenne- 
eigenschappen  uit  nabije  veld  data  besproken  en  onderzoeksgebieden  voor 
de  nabije  toekomst  worden  kort  genoemd. 
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INTRODUCTION 


Since  the  1970s  there  has  been  an  increasing  interest  in  near- field 
measurement  techniques  for  predicting  far-field  antenna  properties.  This 
is  due  to  the  advantages  these  techniques  offer  in  comparison  with 
conventional  techniques  (far-field,  compact  range) [1,  p.xiv;  2,  pp.594- 
596]  . 

In  these  near- field  measurement  techniques,  the  field  of  the  Antenna 
Under  Test  (AUT)  is  detected  in  phase  and  amplitude  by  a  probe  antenna, 
scanning  a  surface  which  is  often  be  only  a  few  wavelengths  away  from 
parts  of  the  antenna  structure.  The  near-field  measurement  technique 
perhaps  most  fully  developed  and  most  easy  to  implement  is  the  planar 
one  and  will  be  discussed  in  the  remainder  of  this  report. 

The  planar  near-field  measurement  technique,  that  includes  the  signal 
reception  in  phase  and  amplitude  by  the  probe,  will  be  described  for  two 
approaches . 

^  In  the  first  approach,  Kerns  (1]  uses  a  scattering-matrix  formulation 
known  from  microwave  circuit  theory.  The  matrices  relate  amplitudes  and 
phases  of  waveguide  modes  and  expansion  coefficients  by  linear  matrix 
transformations . 

In  the  second  approach,  Paris,  Leach  &  Joy  [3]  use  a  Lorentz  reciprocity 
theorem  formulation.  The  AUT  and  the  probe  are  enclosed  in  proper 
surfaces  so  that  a  source-free  volume  is  created.  Then,  using  the 
Lorentz  reciprocity  theorem  for  a  source-free  volume,  an  expression  for 
the  received  signal  is  obtained. 

After  the  description  of  these  two  approaches,  their  equivalence  will  be 
proven  and  the  derivation  of  far-field  antenna  characteristics,  using 
the  results  of  near-field  measurements,  will  be  treated.  _ 
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2  SCATTERING  MATRIX  FORMULATION  APPROACH 

For  the  scattering-matrix  formulation  approach,  the  field  of  the  AUT 
will  be  described  as  an  expansion  of  plane  waves.  Before  this  expansion 
is  given,  the  coordinate  system  of  AUT  and  probe  must  be  given,  as  well 
as  some  wave  parameters  in  these  coordinate  systems. 

2.1  Geometrical  configuration 

Figure  1  shows  the  transmission  system  set-up  for  the  planar  near-field 
measurement  technique  [1,  p.598] : 


receiving 


Fig.l  Transmission  system  set-up  for  planar  near- field  measurement 


Normally  the  AUT  is  operated  as  a  transmitting  antenns  and  the  probe  as 
a  receiving  antenna.  The  aperture  of  the  probe  is  moved  over  a  planar 
surface,  called  the  scan  plane.  Samples  of  the  electric  field,  generated 
by  the  AUT,  are  taken  in  equidistant  points,  that  are  arranged  in  a 
rectangular  grid.  The  properties  of  the  AUT  are  described  in  a 
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rectangular  xyz  coordinate  system  with  its  origin  at  0,  as  shown  in 
figure  1. 

The  scan  plane  is  placed  parallel  to  the  xy-plane  at  a  distance  zo>0- 
The  position  of  the  probe  is  characterised  by  the  point  O' ,  which  has 
the  position  vector  Eg  with  coordinates  (xg,  yg,  Zg).  The  properties  of 
the  probe  are  described  in  a  rectangular  x'y'z'  coordinate  system  with 
its  origin  at  0'  and  its  x’-,  y'-  and  z'-axis  parallel  to  the  x-,  y-  and 
z-axis  respectively.  When  the  probe  is  moved,  its  orientation  and  Zq  are 
kept  constant  while  Xg  and  yg  are  varied. 

The  antennas  are  considered  as  two -port  transducers.  For  the 
transmitting  antenna  one  port  is  placed  in  the  feed  line  at  Sg  while  the 
other  port  is  chosen  to  be  at  the  antenna - aperture .  The  quantities  ag 
and  bg  are  phasor  wave  amplitudes  for  incident  and  emergent  travelling 
waves  of  a  single  waveguide  mode  at  Sg.  The  quantities  a*  and  b*  are 
spectrum  density  functions  for  incoming  and  outgoing  waves  as  defined  in 
the  following.  As  indicated  in  figure  1,  primes  are  used  to  associate 
symbols  with  the  probe. 

2.2  Expansion  in  planar  waves 

The  electromagnetic  fields  in  a  source -free,  linear,  homogeneous, 
isotropic  medium  must  satisfy  the  Maxwell  equations  (an  exp(-jwt)  time- 
dependence  is  assumed  1)  ) 


Vx{- 

(la) 

V  X  a  - 

-jw<£ 

(lb) 

V.fi  -  0 

(lc) 

v.fl  -  o 

(Id) 

^  This  is  the  opposite  of  most  electrical  engineering  formulations,  but 
is  used  here  because  it  is  consistent  with  the  original  derivation  of 
the  plane  wave  theory  [2]  and  is  adopted  by  Newell  [4] ,  upon  whose 
papers  much  of  the  work  to  come  will  rely 
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so  that  a  general  plane  wave  will  have  the  form: 


£  -  £(k)exp(+jfe.l) 

H  -  x  A(k)exp(+Jk-l) 


(2a) 

(2b) 


with: 


k2  -  k-k  -  kx2  +  ky2  +  kz2  -  u2m« 
k-A(k)  -  k^  +  kyAy  +  kzAz  -  0 


(3a) 

(3b) 


in  which  equation  (3b)  originates  from  (lc). 

The  propagation  vector  k  is  regarded  as  a  function  of  its  transverse 

components  lc  ,  k  (which  are  chosen  real).  The  z-component  is  thus: 
x  y 


kz  -  ±7; 


( (^v-k2) 
l  -J(V-2+V-k2) 


if  ^Z+ky2  5  k2 

otherwise 


(4a) 


KV+V 

The  transverse  part  of  k  will  be  denoted  jC,  so  that  K  -  kj^  +  kySy  and: 


7  -  (k‘  -  K  ) 


2  a/2 


(4b) 


7  will  be  taken  positive  for  K<k,  negative  imaginary  for  K>k.  In  (4a), 
+7  is  associated  with  a  plane  wave  travelling  into  the  positive  z- 
direction,  -7  is  associated  with  a  plane  wave  travelling  into  the 
negative  z-direction  (of  course  under  the  condition  that  7  is  real) . 

In  virtue  of  the  relation  k-A(k)  “0.  (2)  yields  just  two  linearly 
Independent  fields,  hence  just  two  basis  fields,  for  any  given  k-  The 
basis  vectors  A  (that  also  indicate  the  polarization  of  the  fi-field 
according  to  equation  (2a)  )  are  chosen  to  be  vectors  parallel  or 
perpendicular  to  the  plane  of  k  and  lz ,  which  is  the  plane  of  incidence 
for  a  wave  Incident  on  any  plane  z-constant  (see  figure  2). 
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Fig. 2  Plane  of  incidence  for  a  wave  Incident  on  plane  z-0 


The  basis  fields  A  are  set  up  with  the  transverse  unit  vectors: 

-  K/K  (5a) 

K2  -  *  *1  (5b) 

which  are  respectively  in  and  perpendicular  to  the  plane  of  Jc  and  gz. 
Both  vectors  and  &2  are  *-n  the  xy-plane.  They  are  illustrated  in 
figure  3. 
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Fig. 3  Propagation  vector  and  transverse  unit  vectors 

The  basis  vectors,  that  must  satisfy  *-  are  chosen  to  be: 

T  K7'^4z,  which  is  parallel  to  the  plane  of  incidence,  and  A^  -  £2 > 
which  is  perpendicular  to  the  plane  of  incidence.  Substituting  and  A^ 
in  (2)  gives: 


El1  -  [&!  *  Ky^tjjJexpC+jk.i;)  (6a) 

Hi1  -  ±Hi*z  x  £iexP(+jk'E) :  Hi  -  “‘/t  (6b) 

&2 ±  -  £2e*P(+Jk-X)  (6c) 

H2±  -  [±h2»z  x  i2  +  *(“»»>  "1*z]«xP(+Jk-i>: 

H2  -  t/(<<v)  (6d) 
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in  which  the  plus-sign  indicates  a  wave  travelling  into  the  positive  z- 
direction  and  the  minus -sign  indicates  a  wave  travelling  into  the 
negative  z-direction.Note  that  (E^~.  represents  a  TM-wave  and  (E2i, 

Ji2~)  represents  a  TE-wave. 

The  transmitted  field  Ed)  of  an  antenna  in  the  region  z>0  for  a 
specified  propagation  vector  fc,  can  now  be  given  as  a  combination  of 
fields  and  E2+’  So  c^e  Plflne  wave  spectrum  expansion  is  given  by: 


E(I) 


00  CO 

J  |  [b1(kx,ky)E1+(K,i)  +  b2(kx,ky)E2+(K,E)]dkxdky 


(7) 


1  2 

in  which  the  coefficients  b  (k  , k  )  and  b  (k  , k  )  are  referred  to  as 

x  y  x  y 

spectrum -density  functions  of  outgoing  waves. 

In  order  to  find  an  expression  for  the  signal  bp' ,  received  by  the  probe 
(see  figure  1) ,  it  is  desirable  to  express  the  radiated  field  E(r)  as 
given  in  equation  (7)  in  terms  of  Ei+(K,E’)  and  E2+(K>E')  associated 
with  the  x'y'z'  coordinate  system.  To  do  this,  use  is  made  of  the 
relation 


£  -  E0  +  I'  (8) 

that  describes  the  relation  between  xyz-  and  x’y'z' -coordinate  system. 
This  relation  is  found  by  inspection  from  figure  4. 


Fig. 4  Derivation  of  translation  theorem 


With  (8)  substituted  in  (7)  it  is  found  that: 

Ei+(K,E)  “  Ei+(lC,£,)exp(+jk.;C0);  i  -  1,  2  (9) 

and  with  (9)  substituted  in  (7),  the  field  incident  upon  the  probe  is 
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00  CO 

E(r')  -  |  J  [b1(kx,ky)e+Jk  K°  E^a.E')  + 

-  CO  .CO 

b2(kx,ky>e+Jkr°  E2+(K,r')Jdkxdky  _ 

00  CO 

-  J  |[a*1(kx,ky)E1+(E,r')  +  a'2(kxfky)22+(kx.ky)]dkxdky 


(10) 


where 

a'i(kx,ky)  -  bi(kx,ky)exp(+jk.£0) ;  1  -  1,  2  (11) 

are  the  spectrum  density  functions  of  the  waves  incoming  on  the  probe. 
2.3  Coupling  equation 

A  two-port  microwave  network,  as  given  in  figure  5,  can  be  described 
with  a  'scattering-matrix',  that  relates  output  wave  quantities  to  input 
wave  quantities  of  the  two-port: 


b0  "  s00a0  +  S01al  (12a) 

bl  “  S10a0  +  Sllal  (12b) 

in  which  Bq,  a^,  bp  and  b^  are  incident  and  emergent  wave  amplitudes  on 
the  two  terminal  surfaces,  0  and  1. 
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Sqq  and  are  the  reflection  coefficients  of  the  two-port;  and  S^q 
are  the  transmission  coefficients. 

The  scattering-matrix  formulation  is  not  restricted  to  two-ports,  it  can 
be  applied  to  any  N-port  microwave -network  (N  -  2,  3,...). 

An  antenna  is  viewed  as  a  multiport  transducer,  having  only  one  input 
port  and  for  each  polarisation  and  direction  in  space  one  output  port. 
This  transducer  transforms  wave  amplitudes  in  a  closed  transmission  line 
system  to  an  angular  spectrum  of  plane  waves  in  the  space  system  and 
vice  versa  (see  figure  6)  [5,  pp.5-6). 


Fig. 6  Antenna  as  multiport  transducer 


The  scattering-matrix  formulation,  relating  input  and  output  quantities 
as  defined  above,  is  given  by  [1,  pp. 602 -603;  5,  p . 7 ] : 
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«  2 


b0  -  r0a0  +  J  J  }  R1OSt.ky)ai(kX'Vdkxdky 


(13a) 


—  1-1 
b^l^.ky)  -  Ti(kx,ky)a0  + 


»  2 


|  f  }  sl,J<>Sc-ky:  1x-VaJ(1x-Vdl: 

>  -®  j-1 


xdly 


(13b) 


with  1  referred  Co  as  polarization  index. 

rQ  corresponds  to  the  reflection  coefficient  at  Sq  looking  toward  the 

antenna  (see  figure  1) ,  R*  and  T*  are  associated  with  the  receiving  and 

transmitting  properties  of  the  antenna,  respectively  and  describes 

the  scattering  properties.  Reception  from  a  certain  direction  (kx.ky) 

depends  on  the  scattering  in  all  directions  (1,1)  from  an  object  in 

x  y 

front  of  the  antenna.  When  primes  are  used,  the  same  relations  are  true 
for  the  probe. 

Now  it  is  possible  to  derive  a  relation  between  the  signal  bQ' (£q) , 
received  by  the  probe,  and  the  input  signal  a0,  delivered  to  the  test 
antenna.  Let  be  the  reflection  coefficient  when  looking  from  the 
probe  into  its  load,  then 


ri'V 


(14) 


To  simplify  the  analysis,  multiple  reflections  between  AUT  and  probe  are 

neglected.  This  means  that  a^(l„,l„)  is  assumed  to  be  zero  and  thus  from 

x  y 

(13b): 


bi(kx,ky)  -  T1(kx,ky)a0  (15) 

Substitution  of  (15)  in  (11)  gives: 

*’l(kx’ky)  *  a0T1<kx,ky)exp(+Jk-I()) 


(16) 
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2 

D(kx-V  "  }  R,i<WTl(kx,ky)  (17d> 

t-1 

DCl^.ky)  is  referred  to  as  the  coupling  product  of  the  receiving  and 
transmitting  properties  of  antennas.  Equation  (17b)  is  known  as  the 
transmission  equation  (4,  p.10}. 
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Equation  (17c)  is  recognized  as  a  Fourier  transform  of  0(1^, ky).  So  the 
coupling  product  can  be  found  by  taking  the  Fourier  inverse  of  (17c). 
Using  the  Fourier  transform  pair  as  stated  by  Papoulis  [6,  p . 1 ] : 


f(t)  -  (2»r)-*  F(u)exp(j<«>t)du>;  F(o>)  -  f (t)exp( - ju>t)dt  (18) 


Taking  the  Fourier  inverse  of  (17c)  yields: 


D(kx.kx) 


1 


*  Vii 

- 7^ - 

4*a  n 


bQ’  (r0)exp(-jjc.i0)dx0dy0 


(19) 


When  the  receiving  characteristic  R'i(kx,ky)  (i  -  1,  2)  is  known, 
Ti(kx,ky)  (i  -  1,  2)  can  be  obtained  with  (19)  and  (17d)  2). 

With  Ti(kx,ky)  (i  -  1,  2)  the  far  field  characteristics  can  be  found  as 
will  be  shown  in  a  forthcoming  chapter.  Now  some  compact  vector- 
formulations  for  equations  derived  in  this  paragraph  will  be  given. 

2.4  Vector  formulation 

The  equations  derived  in  the  preceding  paragraph  will  now  be  written  in 
a  more  compact  vector- formulation  that  corresponds  to  the  formulation 
used  in  the  second  near-field  approach  and  that  facilitates  the 
derivation  of  certain  reciprocity  theorems. 


^  Note  that  in  the  theory  derived  thusfar  only  one  assumption  is  made: 
multiple  reflections  between  AUT  and  probe  are  neglected.  Further, 
due  to  the  presence  of  the  probe  characteristic  R'1  in  the  above 
equations,  it  is  sometimes  stated  that  probe  correction  of  the 
measured  results  is  carried  out. 
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Equation  (7)  can  be  written  as: 


fe(E)exp(+jfc.x)dE 


(20) 


with  fe(£)  following  from  (6)  and  (7): 


fe(E)  -  b1(K)  -  K7'1ft2]  +  b 2(K)JT2  - 

-  b1(E)(k/7]1|(k)  +  b2(K)A  j.(k> 


(21) 


The  unit  vectors  ft|(jj)  and  ft  j.(E)  are  given  by  [2,  p .  95 ]  : 


A|(k)  -  £2  x  ftk:  *k  "  EA  (22a) 

A  ±(k)  -  E2  (22b) 

For  K<k,  where  ft|  and  ft  i  are  real: 

A|(k)-A#(k)  (23a) 

A  J.(k)  -  A*(k)  (23b) 

For  K>k,  ft|  and  ft^  become  complex,  but  remain  unit  vectors  in  the  sense 
A|A|  -  AkAk  -  1  [2,  p.95). 

In  vector  rotation  the  coupling  equation  becomes: 


D(vy  -I'w-xm  (24> 

It  can  be  shown  (appendix  A)  that  for  a  reciprocal  antenna: 


E’(E)-KE)  -  -V?'1 1'(-E).KE);  «?  -  ,/(«/**)  (25a) 

X.'  (E)  •t(E)  -  [l(«/M)1/2/(»»o's)]S’(-E).t(E)  (25b) 
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I(£)  and  £(JC)  are  transverse  vectors ;  they  Include  only  the  components 
transverse  to  the  z-axls.  The  z-component  of  the  complete  transmitting 
and  receiving  vector,  can  be  found  from  the  requirement  that  the 
complete  vector  must  be  orthogonal  to  the  direction  of  propagation  [4, 
p.llj. 


!(£>•£  -0  (26a) 
t(K).£-0  (26b) 

'Complete'  vectors  are  described  by  small  letters,  transverse  vectors  by 
capitals.  The  reason  why  equation  (25)  is  derived  will  become  evident  in 
the  next  chapters . 
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3  RECIPROCITY  THEOREM  FORMULATION  APPROACH 

In  measuring  the  field  of  the  AUT,  the  primary  field  is  altered  by  the 
probe  and  a  true  measure  of  the  original  field  is,  consequently,  never 
attained.  However,  probe  effects  can  in  fact  be  approximated  if  no 
mutual  coupling  between  antennas  is  assumed.  It  is  the  purpose  of  this 
chapter  to  derive  an  expression  for  the  probe  compensation,  by  means  of 
the  well-known  Lorentz  reciprocity  theorem. 

3.1  Geometrical  configuration 

The  transmission  system  set-up  for  the  planar  near-field  measurement 
technique  is  shown  in  figure  7  [3,  p.375]. 
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Fig.  7  Probe  compensation  geometry 


The  field  generated  by  testantenna  A  is  to  be  measured  by  moving  the 
probe  antenna  B  over  the  plane  surface  defined  by  z-Zq . 

Let  the  source-free  volume  V  be  bounded  by  the  closed  surface  2, 
consisting  of  the  infinite  plane  Sg  at  z-a,  the  surface  at  infinity  S^, 
and  the  closed  surface  Sg,  surrounding  the  probe  antenna  B,  operating  at 
the  same  frequency  as  antenna  A,  but  not  necessarily  at  the  same  time. 
£3Sf  Has  represent  the  field  scattered  by  the  test  antenna  A,  when 
generator  B  is  activated;  Egs,  Hgs  represent  the  field  scattered  by  the 
probe  antenna  B,  when  generator  A  is  activated.  The  current  densities 
Ja.  If  ^as '  ^bs  are  deflned  correspondingly. 
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3.2  Lorentz  reciprocity  theorem 

Since  the  volume  V  contains  no  sources,  it  follows  from  the  Lorentz 
reciprocity  theorem  that  [7,  pp. 24-25]: 

[f  KEa+Ebs>x<Hb+Has>  '  (Eb+Eas)x(Ha+Hbs)  1 -n<iS  -  0  (27) 

JJs 

The  integral  vanishes  identically  over  the  surface  at  infinity,  S^.  This 
is  true  because  the  field  is  a  spherical  TEM  wave  on  Sw  [8,  pp.92-93] 
and  the  integrand  in  (27)  can  be  written  as: 


[k/(uM)][(Ea+Ebs>xlnx(V*as>>  '  CEb+Eas)x(Ex(Ea+5bs>,]  »  <28) 

since  for  a  TEM  wave  [9,  p.63]: 

H  -  xn  x  E;  n'1  -  JW O  -  V(^)  (29) 

With  the  vector  identity  &x(£x£)  -  (&-£)£  -  (A-£)£  .  equation  (28) 
becomes : 

[k/(^)lUEa+Ebs).(Eb+Eas))n  -  l (Ea+Ei,s) -n) (Eb+Eas)  + 


-  l(Eb+Eas).(Ea+Ebs))n  +  ( (Eb+Eas) .n)(Ea+Ebs)] -D  (30) 


This  expression  vanishes  identically,  since  E’B  “  0  on  S^,. 
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The  surface  integral  over  Sg  is  best  evaluated  by  introducing  the  unit 
surface  normal  db  -  -n.  If  the  volume  Vg  is  linear  and  isotropic,  the 
Lorentz  reciprocity  theorem  leads  to  [7,  pp. 24-25): 

|f  [<Ia+Wx<VHas>  ‘  (Eb+2as)*(Ha+Hbs)]-(-BB)dS  - 

SB 

-  -f  [  f  [(Et+IW-Ua-^bs)  *  (Ea+£bs>-(4+4s)ldv  <31> 

J  J  JvB 

By  definition,  the  scattered  field  Egs  and  the  currents  and  Jas  are 
zero  troughout  the  volume  Vg.  To  facilitate  the  analysis  of  (31),  an 
approximation  is  made  now:  The  scattered  field  £as  is  neglected  compared 
to  Eb  This  means  that  the  right  side  of  (31)  can  be  expressed  as: 

f  f  f  [Eb-Jbs  *  Ea  ib)^  “  PB<*0>  <32> 

\ 

This  integral  is  proportional  to  the  open-circuit  received  voltage  of 
the  probe  [3,  p.376;  7,  pp. 94-98;  16,  p.1483],  thus  Pg(£g)  represents 
the  measured  signal  within  a  constant  of  proportionality  (see  equation 
(47)  ). 

Equation  (27)  can  now  be  written  as: 

j  J  [ (Ea+Bbs)x(Hb+Has)  -  (Eb+Eas>x(Ha+Hbs>)-ndS  -  PB(lo)  (33) 

J  Jso 
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In  this  integral,  the  terms  involving  products  of  primary  and  scattered 
fields  vanish  identically.  For  example, 

|  }  tEa  x  Has  '  Eas  x  Ha]  -HdS  -  jj  x  Has  -  Ias  x  HJ  .ndS  - 
s0  S0+S« 

-  [  f  [  tEas4  '  Ea.Jas]dV  -  0  (34) 

Jv+vB 

The  last  step  is  true  because,  within  V  and  V„ ,  the  currents  J  and  J„„ 

a  —a  —as 

are  zero  by  definition;  the  first  step  is  permissible  because  the 
contribution  from  the  surface  at  infinity  tends  to  zero. 

Terms  in  (33)  which  involve  only  scattered  fields  will  be  small  compared 
to  terms  only  involving  primary  fields;  hence  they  will  be  neglected.  So 
(33)  becomes: 

J  J (Ea  x  Hb  -  Eb  X  Ha] -ndS  -  PB(r0)  (35) 

s0 

Further  reduction  of  (35)  can  be  effected  by  expanding  the  primary 
fields  in  terms  of  their  wavenumber  spectra  over  the  surface  Sq  (z-a) . 
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With  a  time -dependence  according  to  exp(-jwt): 


Ea(x,y,a)  -  £(E)exp(+jk.r)dE 


Ha(x,y,a)  -  Is  x  £(JC)exp(+jls.E)d2C 


IS-£(K)  -  0 


00  00 

^(K')  -  J  J  g(K'  )exp(+jjc'  )dJC' 


Hb(K')  -  (Ufi) 


00  00 

i  j*'- 


*<r)exp(+jk'.E')d*' 


k'.g(K')  -  0 


and  (35)  becomes: 


00  00  00  00  00 


J  J  J  J  J  t£(K)x{Js'xg(E'))  -  g(E')x{kx£(E))]. 


-00  -  OO  -  CO  -CO  -CO 


(-4z)exp(+jls.r  +  jls' .r')dE'dJ5dxdy  (37) 
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The  evaluation  of  (37)  is  rather  involved  and  will  be  presented  in 
appendix  B  [10].  It  appears  that  (37)  can  be  written  as: 


[4x2/(cd/i)J J  j  [f(I)x(kxg(-K))  +  *(-K)x{]sxf(E))J.4zexp(+jfc.Io)  - 

-  00  -00 

-  PB(lo>  (38) 

Using  the  vector  identity  Ax(l x£)  -  -  (A.B)C  and  equations 

(36c, f)  yields: 


[8  *2/(w/i)] 
-00 


kzf(K).g(-K)exp(+jlc.  £<))<« 


(39) 


Equation  (39)  is  recognized  as  a  two-dimensional  Fourier  transform. 
Taking  the  inverse  over  the  plane  z-z0  gives  (see  equation  (18)  ): 


f(K) •*(-*) 


1 


00  00 


PB(E0)exp(-jk.£0)dx0dy0 


(40) 
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4  EQUIVALENCE  OF  BOTH  APPROACHES 

From  inspection  of  equations  (20)  and  (36)  for  the  scattering-matrix 
formulation  approach  and  the  reciprocity  theorem  formulation  approach, 
respectively: 


E(£) 


fe(E)exp(+jfe.E)dJC 


Scattering-Matrix  Formulation  Approach 


(41a) 


ECt) 


f(K)exp(+jfe.E)d|C 


Reciprocity  Theorem  Formulation  Approach 


(41b) 


it  follows  that 


f(E).*(-K)  -  fe(E).fe'(-E)  (42a) 

because : 


E(l)e*P(-jk'l)dx<iy 


(42b) 


-  f  (K) 
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and  this  substituted  in  (45)  with  k  -  u7(<m)  gives: 


Pb(£0>  -  8»2a02^0 


oo  eo 

»J  Jl'(K)- 


t(I)exp(+jk. ig)djC  - 


o  2 _  2  1  :  ro'ri'  f  ao 

Sir  aQ  r,0  -  ( 

a0  i  1  "  r0  rl 

CO  00 

|  J  t'(K).t(K)exp(+jlc.i:)dK  |  - 

-CO  -00 

-  8ir2a0iJ0(l  -  ro'rj'JbQ'  (£g) 


with  bg' (Eq)  given  in  (17b, c). 

Substitution  of  (44),  (46)  and  (47)  in  (40)  gives: 


(47) 


a 


2 

0 


IQ* 

7  JU/V) 


t(K).r’(K) 


l 

(2s)2 


8ir2a0r,0(l 


W 


). 


b0’(E0)exp(-jk.l0)dx0dy0 


(48) 


or,  with  (24): 


Dd^.ky) 


<WWp)  a  •  iyry) 


k(2*)“ 


Jv 


(l0)exp(-jjE.t0)dxod>r  " 


oo  m 

^  i 


bo'(l0)exp(-jl5.K0)<lx0dy0  (49) 
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with  Dd^.ky)  -  t(K).r'(K)  -  KK).E'(K). 

This  is  the  same  result  as  stated  in  equation  (19)  and,  consequently, 
both  approaches  are  equivalent. 
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5  FAR- FIELD  ANTENNA  CHARACTERISTICS 

Now  that  Is  proven  that  both  approaches  lead  to  the  sane  results,  a 
choice  In  favour  of  one  of  the  two  approaches  has  to  be  made .  Chosen  Is 
for  the  scattering-matrix  formulation  approach.  Maybe  this  approach  is  a 
little  less  transparent  than  the  reciprocity  theorem  approach,  but  it 
seems  that  the  scattering-matrix  formulation  approach  is  the  most  widely 
adopted  one  in  the  area  of  planar  near-field  analysis  [11],  obviously 
due  to  the  pioneering  work  done  and  expertise  available  at  the  U.S. 
National  Bureau  of  Standards  (NBS) . 

In  the  remainder  of  this  chapter,  equations  for  the  far-field,  power- 
gain  function  and  receiving  area  of  an  Antenna  Under  Test  will  be 
derived. 

5.1  Far-field  of  Antenna  Under  Test 

The  starting  point  for  the  derivation  of  the  far-field  formulation  is 
equation  (20): 


E(I) 


fe(K)exp(+jfc-K)dK 


(50) 


With  equation  (43) : 


E(r) 


£(K)exp(+jk.l)dK 


(51) 
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Normalized  to  Che  Input  aignal  (a0-l) ,  the  E-field  becomes: 


E(l) 


KI)exp(+jl5.£)dK 


(52) 


By  the  method  of  the  steepest  descent  (known  also  as  the  saddle  point 
method  of  integration),  or,  alternatively,  the  method  of  stationary 
phase,  (32)  can  be  approximated  in  the  far  field,  that  is  for  large  r. 

As  proven  in  appendix  C  by  the  last  mentioned  method  [12,  pp. 28-38),  the 
far- field  formulation  is  given  by: 


+jkr 

E(l)  =  -J2»k z0£(E)— - ;  k  o  -k  — 

r-~°  r  r 

+jkr 

-  -j2xkcos0£(E)-- - ;  k  0eR  (53) 

r  z 

This  equation  differs  from  the  far- field  given  by  Newell [4,  p.19]  by  a 
factor  2ir,  due  to  an  applied  normalization  that  cancels  out  this  factor 
[2,  pp. 271-274],  The  equation  as  stated  above  is  also  found  in  [17, 
p.ll]  and  will  be  used  in  the  remainder  of  this  report. 

5.2  Power-gain  function 

The  power-gain  function  is  given  by  [8,  p.43]: 


Gain 


radiation  Intensity  U(0,^) 

4* -  -  4* - 

total  input  power 


(54a) 
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in  which  Che  radiation  Intensity  is  given  by  [8,  p . 28 ] : 


T*" 

U(«,*) - —  U  j£(i) |2;  r,  - 


(54b) 


rj  is  the  characteristic  admittance  of  the  medium. 
The  total  input  power  P^n  is  given  by  [13,  p.536] : 


Pin  -  1/2  Re(V0I0  ) 
With  [14,  p. 24] : 


(55) 


v0  -  a0  +  bQ 
Iq  ”  'I o(aO  '  b0^ 


(56a) 

(56b) 


ij0  is  the  characteristic  admittance  of  the  antenna.  So  the  input  power 
is: 


Pin  “  1/2  i70Re((a0  +  b0)(a0 


bn*))  - 


“  1/2  ijqRc  { &Q&Q  *  b0bQ  +  3q  bQ  -  «QbQ  )  — 

-  1/2  ijg(  I a0 1 2  -  I b0 1 2 )  (57) 

off  ^ 

since  aQ  bQ  -  agbQ  is  pure  imaginary. 

With  bg  -  Tgag  (Tg  is  the  reflection  coefficient  when  looking  into  the 
antenna) : 


pin  “  V2  *»0l“ol2<1  *  lrol2>  (58> 

The  power  gain  is  found  by  substituting  (53),  (54b)  and  (58)  in  (54a) 
with  aQ  -  1  (because  of  the  normalization  of  equation  (51)  ). 
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For  the  gain  is  found  (see  also  {1,  p.606]): 

o  9  I £(E)  |  ^*7  , 

G(E)  -  16*3k2cos20— - - j - - ;  v  -  7(*/m)  (59) 

-  lror>*»o 

This  equation  differs  from  that  given  by  Newell[4,  p.19],  due  to  the  2n 
difference  in  the  far- field  equation. 

5.3  Receiving  effective  area 

The  receiving  effective  area  of  an  antenna,  that  is  smaller  than  the 
physical  aperture,  is  related  to  the  power-gain  according  to  [15, 
pp. 178-179]: 


«(K)  -  [A2/(4*)]G(g);  X  -  2*A  (60) 

To  express  the  receiving  effective  area  totally  in  receiving 
characteristics,  use  is  made  of  equation  (46),  and  together  with  (59) 
and  (60)  this  yields: 


o(K) 


A216ir3k2cos20  |e'  (-K)|2q02k2q 

4*<i- 1  r0 1  2  >i70T'^’?2 


16*4- 


ll'(-K)l2>70  . 

u  -  lr0l  in  ’ 


n  -  A«/p) 


(61) 


The  difference  between  the  equation  (61)  and  that  stated  by  Newell  [4, 
p.19]  is  caused  by  the  2w-difference  in  the  far-field  equation. 
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6  RESEARCH  ITEMS 

Now,  that  the  theory  of  near-field  measurement  is  understood,  attention 

can  be  paid  to  the  following  items: 

co  -  and  crossoolar  reception:  The  elctric  field  vectors  of  AUT 
and  probe  can  be  split  in  two  orthogonal  components .  In  one 
measurement  the  probe  will  couple  primarily  to  one  of  the 
components,  called  the  copolar  component,  the  other  component  is 
referred  to  as  the  crosspolar  component.  In  a  second  measurement 
with  the  probe  90  degrees  rotated,  it  will  couple  primarily  to 
the  other  component,  then  referred  to  as  copolar  component.  In 
order  to  obtain  the  desired  components  in  two  measurements,  first 
an  appropriate  definition  of  crospolarization  has  to  be  stated; 
samu-it.E  algorithm:  First  the  distance  between  scanning  plane  and 
,Trr\e  has  be  chosen.  Then  the  sample  spacing  and  the  number  of 
sample  points  must  be  selected  to  include  the  details  of  the 
field  in  the  aperture; 

fflea?uxl.ng -flPtiVg  Ohtennas_in  transmitting  and  receiving  mode;  In 
the  theory  derived  before  it  was  assumed  that  the  AUT  is 
operating  as  transmitting  antenna  and  the  probe  as  receiving 
antenna.  For  measuring  an  antenna  in  the  receiving  mode,  these 
roles  have  to  change.  It  is  very  likely  that  the  theory,  with 
some  little  modifications,  will  hold  true  for  this  situation; 
error  analysis:  In  order  to  draw  conclusions  from  the  results 
obtained  from  near- field  measurements,  a  thorough  understanding 
of  possible  error  sources  and  their  influence  on  the  total 
results  must  be  available .  Of  the  three  possible  error 
categories:  theory,  numerical  calculations,  and  measurement,  the 
latter  will  be  responsible  for  most  of  the  errors  involved; 
pulsed  measurements:  The  consequences  of  pulsed  measurements  for 
the  near- field  theory,  that  is  derived  under  the  assumption  of 
unmodulated  Continuous  Vave  operation,  have  to  be  investigated. 
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7  CONCLUSIONS 

A  successful  near-field  antenna  measurements  program  requires,  besides 
careful  measurements  and  computer  processing  of  the  data,  a  reasonable 
understanding  of  the  theory  involved. 

This  theory  has  been  covered  in  detail  in  this  report,  following  two 
approaches,  and  can  be  used  as  a  guide  in  the  understanding  of  the 
theory  of  planar  near-field  measurement. 


Ir.  J.G.  v.  Hezewijk 
(Project  Leader) 


Ir.  H.J.  Visser 
(Author) 
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APPENDIX  A  RECIPROCITY 


To  prove  is: 


'?0k£'(K)  *  1  “  7(«/m)  (A.l) 

In  order  to  show  the  validity  of  (A.l),  first  the  following  so-called 
reciprocity  lemma  is  proven: 

r,0(a0'b0  -  a0b0')  -  (a'(E)b(K)  -  a(£)b'  (JC)  ]q(K)dK  (A.  2) 

JK 


A.l  Reciprocity  lemma  (relationship  between  transverse  receiving 

and  transmitting  characteristics) 

Equation  (A. 2)  relates  the  electromagnetic  fields  on  the  AUT-aperture  to 
the  electromagnetic  fields  on  the  probe -aperture .  To  prove  (A. 2),  first 
the  field  description  on  an  arbitrary  aperture  S,  perpendicular  to  the 
z-axis  (see  figure  1),  is  given. 

A. 1.1  Transverse  field 

The  electric  and  magnetic  fields  Et  and  Ht  must  satisfy  the  following 
Maxwell  equations  on  S  (an  exp(-jwt)  time-dependence  is  assumed): 


J«^HC 

v  X  Ht  -  - jw«Et 


(A. 3a) 
(A. 3b) 
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In  the  rectangular  coordinate  system  xyz  of  figure  1,  the  7-operator  can 
be  divided  into  a  transverse  and  an  axial  part  according  to  [Al,  p . 67 ] : 

“  &x  S/Sn  +  Ay  S/Sy  (A. 4a) 

7z  "  &z  S/Sz  (A. 4b) 

With  (A. 4)  and  the  fact  that  £t  and  Ht  do  not  possess  a  z-component, 
(A.l)  can  be  written  as: 


ft2  x  S/S z  £t  -  jupHt  (A. 5a) 
a2  x  S/Sz  Ht  -  -jw<£t  (A. 5b) 
Vt  x  £t  -  2  (A. 5c) 
Vt  x  gt  -  2  (A.5d) 


The  last  two  equations  indicate  that  £t  and  Ut  can  be  written  as, 
respectively  [Al,  p.68]: 


Et  -  g1(z)7t«  (A. 6a) 

Bt  -  g2(z)7t#  (A. 6b) 

Equations  (A.3a,b)  implicate: 

S2/Sz2  £c  +  k2£t  -  2;  k2  -  u2tp  (A. 7) 

and  with  (A. 6a)  substituted  in  (A. 7),  it  is  found  that  g^(z)  can  be 
written  as : 


gj^z)  -  A+exp(+jkz) 


(A. 8) 
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So  for  the  electric  field  is  found: 

Et  “  A±7t#(x,y)exp(±jkz)  (A. 9a) 

and  with  (A. 5a): 

Ht  -  *V&Z  x  A±Vt*exp(±jkz) ;  rj  -  J(t/n)  (A.  9b) 

A. 1.2  Two-port  description  of  AUT-probe  system 

The  Aut-probe  system  is  modelled  as  a  two  port,  see  figure  A1 
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Fig.  Al  Two-port  model  of  Aut-probe  system 


Incident  and  reflected  voltage  waves  at  port  n  (n-1,2)  are  given,  in 
amplitude  and  phase,  by  VR+  and  Vn‘ ,  respectively.  Incident  and 
reflected  waves  are  normalized  according  to  Kerns  &  Beatty  [A2,  p.2A] 
Al) 


-  V, 


+  . 


b  -  V* 
n  n 


(A. 10a) 


Al) 


The  normalization  is  somewhat  arbitrary.  The  results  are  normali¬ 
zation-independent,  as  long  as  the  normalization  is  such  that 
I/V  -  t)  for  an  incident  or  reflected  wave  [A. 5,  pp. 525-526). 
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Voltage  and  current  on  the  reference  planes  are  given  by: 


V  -  V*  +  V  '  -  (a  +  b„)  (A. 10b) 

nnn'nn 

I„  -  C„(v  +  -  V  •)  -  C_(a_  -  b  )  (A. 10c) 

n  n'  n  n  '  n  n  n 

Vn  and  In  are  introduced  to  facilitate  the  calculations;  they  are  not 
available  for  measurement. 

The  characteristic  admittance  of  the  ntb  port  is  connected  in  a  simple 
way  with  Cn  [A. 6,  pp. 146-147],  To  show  this  connection,  let  bn-0 ;  then: 


V 
I 

V 


n 

n 

n 


Vn 

^n 


C 


n 


When  the  antennas  are  matched  to  free -space: 


(A. 11a) 
(A. lib) 
(A. 11c) 


Cn  "  ’'n  “  •»  (A.  lid) 

The  above  is  now  applied  to  the  theory  in  paragraph  A. 1.1  by  relating 
the  voltage  to  the  electric  field  and  the  current  to  the  magnetic  field: 


(A. 12a) 
(A. 12b) 

Substitution  of  (A. 10)  and  (A. 11)  in  (A. 12)  gives: 


fit  -  w 
fit  -  W 


fitn  -  <an  +  W  -  <“n  +  V*n 
Htn  -  <an  -  “  <*n  '  bn^ 


(A. 13a) 
(A. 13b) 
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(A. 9)  yields: 

Btn  -  •»**  x  *tn  (A14) 

or,  with  (A. 13): 


a_  +  b_ 
n  n 


(A. 15) 


Substitution  of  (A. 10)  and  (A. lid)  in  (A. 15)  gives: 


V 

K  -  1  ~T~7T~  x  fin 

W 

o  V„ 

-  >72-rL-ft 2  x  4n  - 

in 

-  iV1!*  x  ” 

"  9fiz  x  fin 


(A. 16) 
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A. 1.3  Reciprocity  theorem 

Consider  the  two-antenna  system  of  figure  A2  [A3,  p.95] 


Fig.  A2  Two -antenna  system 


In  the  source -free  region,  bounded  by  the  surfaces  S,  Sq,  Sc,  Sc'  and 
the  following  Lorentz  reciprocity  theorem  is  valid: 

(1*1'  -  E'  x  H)  -nds  -  o  (a.  17) 

X  consists  of  the  (perfect)  conducting  surfaces  Sc  and  Sc',  the 
radiating  surfaces  S  and  Sq  and  the  surface  at  infinity  Sw. 
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The  contributions  from  Sc  and  Sc'  to  (A. 17)  are  zero,  because  Et  -  2  on 
Sc  and  Sc'.  On  SB  is  H  -  qAr  x  £  .  end.  consequently: 


j  |  (ExH'  -  E'xH).ndS  -  qj  |  (ExArxE-  -  E'x&rxE).firdS  - 
-  qj  |  f*r(E-E')  -  E' (E-ftr)  -  *r(E'-E)  +  £(£' -£r) ] • 


ardS 


-  0 


So  for  the  system  of  figure  A2,  (A. 17)  results  in: 


(A. 18) 


Et'xflt).-AzdS 


-  J  J  (Et '  xHt 


EtxHt’)-fizdS  (A. 19) 


The  primes  in  this  equation  are  associated  with  the  Probe  (see  figure 
A3). 
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On  the  antenna  apertures  S  and  Sq,  the  fields  Et  and  only  depend  on 
the  transverse  part  of  (E)  and  the  transverse  coordinates  (x,y).  The 
integrand  on  (A. 19)  can  be  written  as  (with  reference  to  figure  A3): 


Etn*Htn'  -  Etn'xHtn  -  f  {[»n(E)+bn(K) )  [«„'  (-»  -b„'  (-» ]  (4^'  )  + 

E 

-  [«n'(-K)  +  bn'(-E))[an(E)  -  bn(E)](en'xhn)|dK 
|>j(E)|[an(K)an'  (-g)-bn(E)bn'  (-g)+an'  (-JC)bn(K)-an(K)bn-  (-K)  ]  . 


(A. 15) 


E 

.(Snx4zx4n,)  + 

-  [an(E)an' (-K)-bn(E)bn' (-K)-an’ (-K)bn(E)+an(E)bn' (-K) ] . 


(Vx*zx*n> 


jdE  - 


-  |  <?(E)|[an(K)an'(-E)-bn(E)bn'(-E)+an'(-E)bn(K)-an(K)bn'(-K)]. 


(«z<V*n'>  -  4n'<Sn-4z>>  + 

-  [an(E)an' (-E)-bn(E)bn* (-E)+an(E)bn’ (-E) (-E)bn(E) 1 ■ 
•t*z(V  *n>  -  *n<V  “ 

-  2ftz[  fKE)(an'(-K)b(E)-an(E)bn’(-E>HvV]dE  <A  20) 


The  dot  product  t*n-*n  ']  in  the  integrand  of  (A. 20)  Illustrates  the 
Influence  of  the  polarization. 
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Returning  to  (A. 19),  taking  reference  plane  Sq  In  the  waveguide  behind 
the  probe,  where  K  “  2  (see  figure  A3)  and  taking  reference  plane  S  on 
the  ATJT  aperture,  equation  (A.  19)  becomes: 


,*o[ao'Vaobo']tso-V!dxodyo  ” 
Js0 


[a'(-E)b(£)  + 

K 


-  a(E)b *(-K)  ]•>(*)] 


} 


[ft.e' ]dxdy 


(A. 21) 


or: 


90(a0'b0-a0b0’) 


I  I*'* 


’  dxdy 


40-fiO,dxOdyo 

Js0 


ta'(-E)b(E)-a(K)b'(-K))^(K)dK 

S 


(A. 22) 


Including  both  polarizations: 


r>o^ao’bo'aobo'  ^ 


ft. ft' dxdy 
JS 

I  4o.4o'dx0dy0 
S0 


I  [ft '  (  -B)  -fe(E)  -ft(E)  -ft’  ( -K)  ]  r,  (E)dK 

E 


(A. 23) 


(Note  the  dot  product) . 

When  the  medium  is  isotropic,  ij  is  independent  of  K  and  can  be  placed 
outside  the  integral. 
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Equation  (A. 23)  is  also  valid  for  the  situation  wherein  AUT  and  probe 
are  interchanged,  with  reference  planes  on  the  probe  aperture  and  in  the 
waveguide  behind  the  AUT. 

A. 1.4  Antenna  reciprocity 

To  derive  an  equation  relating  the  transverse  receiving  characteristic 
&(&)  to  the  transverse  transmission  characteristic  I(K) ,  use  is  made  of 
equations  (13a,b). 

Scattering  equations  for  AUT  and  probe  are  given  by: 


00  QO  2 


b0  -  r0a0  +  J  |  \  Ri(K)ai(K)dK 

-  CD  -  00  t- 1 

00  CO  2 

k^OS)  -  Ti(S)aQ  +  |  |  ^  Si-J(K,I,)aJa()dI( 


-co  J=1 
•  2 


b0'  -  r, 


o'V ’  +  |  |  ^  R'1®*'1®® 

-00  -  00  i— 1 

ao  oo  2 

b'1^)  -  T'l(K)a0'  +  |  |  ^  S'l'J(K,L)a'Ja,)dl 


(A. 24a) 


(A. 24b) 


(A. 24c) 


(A. 24d) 


With  a  proper  choice  of  excitation  the  wanted  equation  can  be  derived. 
Consider  the  excitation: 


a0  -  1;  ■  *(*)  -  0 

■0'  -  0;  a'J(K)  -  S^SCE-l.) 


(A. 25) 


This  means:  The  AUT  is  transmitting  and  the  probe  is  receiving  from  the 
direction  of  the  AUT,  see  also  figure  6,  page  11. 
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With  (A. 22)  is  found: 


J 

£.£'dxdy 

JS 

"JoV  "  •»  | 

I  £0  .£<)'  dXQdyg 
3s0 

I 

£.£'dxdy 

JS 

T 

j  fio-fio'dXo^O 

Js0 

I 

£.£’dxdy 

Js 

'T 

1  ftg.SLo'dxodyo 
Jso 

J 

£.£'dxdy 

Js 

‘7 

1  fig.fig'dXgdyg 
*0 

since : 

'  {  0 

for  i-j 

else 

«<-L-E)  -  |  j 

L  for  Ir-K 

)  else 

ta’Jc-Db^Dldl 


Jv 


b^-K) 


(A. 26a) 


(A. 26b) 


(A. 26c) 
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With  the  excitation  according  to  (A. 25),  (A. 24)  yields: 


b1(K)  “  T1^)  (A. 27a) 

V-R’^E)  (A. 27b) 


so  (A. 26a)  transforms  to: 


|  |  fi.fi'dxdy 

-  <!  f  -  Ti(-E)  (A. 28) 

I  I  fio'£o'dxOdyo 

S0 

When  both  antennas  are  identical,  the  wanted  equation  is  found  to  be: 


<70E(K)  ~  (A- 29) 

since  S  and  Sq  are  the  same  in  magnitude,  see  figure  A4.  It  is  obvious 
that  S  -  Sq  when  the  probe  is  a  waveguide.  For  a  horn  antenna  (see 
dashed  line  In  figure  A4) ,  the  fields  at  the  aperture  can  be  found  by 
treating  the  horn  as  a  radial  waveguide  (dotted  line  in  figure  A4) 
radiating  the  dominant  mode  of  the  waveguide- feed  [A7,  pp. 533-534;  A8 , 
pp.349],  so  also  in  this  case  S  -  Sq. 
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WAVEGUIDE  * 


■>  FREE  SPACE 


& 


^2 


s.=  s 


Fig.AA  Antenna  aperture 


A. 2  Relationship  between  complete  transmitting  and  receiving 

characteristic 

In  analogy  with  (21)  (that  is  related  to  a  transmitted  wave),  the 
complete  transmitting  characteristic  is  given  by  [A9,  p . 1 2 ] : 

id)  -  [•c/7)T1(K)fi  |  (fc)  +  T 2(g)ft  x(k)  (A.  30) 

In  analogy  with  (7),  that  gives  the  transmitted  field  of  an  antenna,  the 
received  field  is  given  by: 


(A. 31) 
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With  (6a)  and  (6c)  the  Integrand  can  be  written  as: 

a1(K)E1'(K,r)  +  a2(K)E(E,r>  - 

-  +  k7_1A2)  +  a2(E)£2}exp(+j]t.£)  (A.  32) 

This  means  that  the  received  transverse  field  Et(£)  is  given  by  [2, 
P-91] : 


Et(r> 


fa1(K)^1  +  a2(g)^2]exp(+jk.r)dkxdky 


A(K)exp(+jk-E)dkxdky 


(A. 33a) 


with: 


A (K)  -  a1(E)E1  +  a2(K)^2 


(A. 33b) 


Now  assume  that  reception  from  only  one  direction  takes  place,  then: 


Et(E)  -  A(E)exp(+jls.£) 


(A.  34) 
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With  (13a): 


m  te  r-2 


o  -  roao +  }  J14Rl<*)*i<Bdk*dky 


CO  CO 

-  r0a0  +  J  |  E(K)  .ACKJdl^dky 


-  r0a0  +  E(K)  -A(IC) 

reception  from  one  direction 

with  A(K)  as  defined  above  and: 


(A. 35) 


£(E)  -  R1^!  +  R2(E)K2  (A. 36) 

To  avoid  confusion  with  vector  A  as  used  in  paragraph  2.2,  A(K)  in 
(A. 35)  will  be  rewritten  £(K) : 


b0  "  r0c0  +  KK)-fi(S)  (A.  37a) 

£(K)  -  c1(E)X1  +  c2(K)g2  (A. 37b) 

Vector  £  is  given  by  (compare  with  (21)  ): 

£(K)  -  +  Ky ”  1az ]  +  c2(E)K2;  kz  -  -7  (A. 38) 

kz  -  -7  because  the  wave  is  incident  upon  the  antenna. 
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With  (5)  and  (22): 


4|(k>  -  [lAK](4zxg)xj£  - 

-  -[l/kK]fex(ftzx£)  - 

-  -[1/kK]  (az(k  K>  -  E(k.4z>)  - 

-  -I1AK](K2AZ  +  tK)  - 

-  -{[KA]*Z  +  [7AHE/K])  - 

-  -[7A](Ex  +  K-r'1^)  (A. 39) 

and  thus  equation  (A. 38)  becomes: 

c(K)  -  - [k/Tf]cX<K)*|  (Is)  +  c2(E)a  jL(k)  (A. 40) 


The  complete  receiving  vector  l(K)  is  found  from  the  requirement  [A4, 
P-125] : 


R<E)-£(S)  -  i(E).fi(E) 


(A. 41) 


(which  states  that  z-components  are  redundant  in  describing  the 
electromagnetic  fields).  Equation  (A. 37)  gives: 


X(E)  -  -ItAIR^E^iOs)  +  R2(E)4  i(E) 


(A. 42) 


This  can  be  checked  by  substituting  (A. 42)  in  (A. 37): 

1(E). fi(E)  -  (  -(7A)R1(E)*|(k)  +  RZ(E)4  j.(k>  )• 

.(  -(R/7)c1(E)4i(k)  +  c2(E)4  x(k)  )  - 

-  R1(E)c1(E)  +  R^(E)c2(E)  - 

-  (  R1(E)E1  +  R2(E)S2  )•(  cl(E)Ei  +  c2(K)E2)  - 

-  R(E)-£(E) 


(A. 43) 
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Note  that  the  following  equation  is  valid  too: 


E(E)KE)  -  t(E).KE) 


(A. 44) 


With  use  of  (A.  29),  the  desired  relation  between  X(K)  and  £(jC)  will  be 
derived. 

According  to  (A. 42): 


r>0ki( K)  -  +  i70kR2(E)a  i(k)  - 

(A. 29)  i  1  i 

-  '♦■^0'yr7Xr7o” lT± ^ (k)  *  10*1210  l(ls)  “ 

-  Vlili/k]  [k/7)T1(-K)a^(k)  -  t,2kI2(-S)«  i(k)  - 

-  -r(^/Y][7A][k/7]T1(K)t|(k)  -  [T/M]kT2®a  i(k)  - 

-  [«e7A][k/7]T1(E)fi|(k)  -  [Tk/(w/0  ]T2(E)ft  j.(k)  - 

-  7*?[k/'r]T1(K)ft|(k)  -  79T2(K)A  l(lt);  r,  -  7(£/m)  (A. 45a) 


In  the  first  step  above,  use  is  made  of: 


-  ui/t,  fi2  -  l/(wn)  (A. 45b) 

Now,  first  the  unity  vectors  A|(k)  end  ft  j_(k)  are  considered. 

With  (22)  and  (5): 

ft  ±(k)  -  E2  -  ftz  x  ^  -  ft2  x  EA 

»  ft  l('k)  -  ftj,  x  -EA  -  -ft  l(k)  (A. 46a) 

ft|(k)  -  E2  x  ftk  -  ft  i.(k)  *  kA 
»  ft| (-k)  -  t  i(-k)  *  -kA  - 
-  -ft  i(k)  x  -kA  -  ft|  (k) 


(A. 46b) 
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Equation  (A. 46)  substituted  in  (A. 45)  yields: 


l0kl(K)  -  -r*7(  [k/7]T1(-K)*|(-k)  +  T2(-K)A  jL(-k)  )  - 

“  7*7t(-K);  “  7<£/m)  (A. 47) 


which  is  the  equation  searched  for. 
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APPENDIX  B  INTEGRAL  EVALUATION  [Bl] 

To  be  evaluated  Is  the  Integral  equation  (37): 


V*o> 


CO  CO  00  00  00  CD 

<*>1  inn 


[£(E)x(k'xa(K')>  -  B(K'  )x(Jcx£(J()  )  J  . 


•  <-42>exp(+jlt.i  +jk'  •  I  ’  )  dKdjC  ’  dxdy  ’  (B.l) 


B.l  Evaluation 

The  vectors  k,  k' ,  E  and  e'  can  be  written  as: 


k  -  kxBx  +  ky*y  +  kztz  (B. 2a) 

k'  -  V*x'  +  W  +  V*z'  (B.2b) 

I  -  Xftx  +  y*y  +  z*z  (B.2c) 

E'  x'»x’  +  y'ty'  +  z'*z'  (B . 2d) 


The  two  rectangular  coordinate  systems  are  shown  in  figure  Bl  (see  also 
figure  7) . 
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and: 


Xq  -  £  -  I'  -  (x+x')4x  +  (y-y')ly  +  (a+b)ftz  - 

'  rOx*x  +  r0y*y  +  r0z*z  <B5a> 

*»  x'  "  r0x'x:  y'  -  yr0y:  b  "  r0z’a  (B  5b) 

The  variables  x  and  y  appear  only  in  the  exponential  term  of  the 
integrand.  Integration  with  respect  to  them  leads  to: 


exp(+jfc-E  +jfc’ .£* )dxdy 


exp  <+ j  kxx+ j kyy+ j kzz+j  kx ’ x '  + 


+Jky'y'+jkz’z'  )dxdy 


-  expl-jk'.^Q  +  ja(kz-kz'))|  J  exp( jx(kx-kx' )  +  jy(ky+ky ’ )dxdy 

(B.6) 

This  is  recognized  as  a  double  inverse  Fourier  integral.  A  single 
inverse  Fourier  integral  looks  like  [B2,p.l9]: 


(1/(2*)] 


I 


exp(+Jyf)dy  -  5(f) 


(B. 7a) 


with  S( f)  the  Kronecker  Delta  function. 
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With  (B. 7a) : 


1 

CO 


expl+jx^-l^' )  )dx  -  2*4(kx-kx’) 


<B.7b) 


Using  this  in  equation  (B.6)  gives: 


expCjfc.i  +  j]c'.£')dxdy  -  (2x)Zexp( -jjt'  .Iq  +  ja(kz-kz'). 


•4(W>*(ky+ky’>  (B.8) 


and  this  result  is  substituted  in  (B.l)  and  integrated  with  respect  to 
kx'  and  ky' .  In  accordance  with  the  sampling  properties  of  the  Dirac 
function  (4),  this  leads  to  a  simple  evaluation  of  the  rest  of  the 
integral  at  kx'  -  1^;  ky'  -  -ky.  It  will  be  noticed  that  the  third 
component  of  fe'  is: 


V  -  -  V2  •  ky'2)  -  ■  K2  *  ky2>  -  kt 

and  thus : 


CB.9) 


k'  -  -k 

Equation  (B.l)  can  now  be  written  as: 


(B.10) 


[4w2/(«fi)] 


iKK)xlkxg(-K) 


+ 


+  g(-K)x{Jsxf(I))J  ■lzexp(+JJj.£0)dX 


(B. 11) 
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APPENDIX  C  METHOD  OF  STATIONARY  PHASE 


Consider  che  function  I: 


I 


f(x,y)exp! jkg(x,y) )dxdy 


(C.l) 


For  a  large  k,  the  phase  of  I  will  vary  rapidly  and,  under  the  condition 
that  f(x,y)  is  a  slowly  varying  function,  I  will  tend  toward  zero  when  k 
approximates  infinity. 

If,  however,  there  are  points  where  g(x,y)  is  stationary,  that  is 
Vg(x,y)  -  0,  the  term  exp{ jkg(x.y) )  will  not  vary  rapidly  in  the 
vicinity  of  these  points  and  the  main  contribution  to  I  will  be 
delivered  in  these  points.  By  expanding  the  functions  in  I  around  these 
points  of  stationary  phase,  an  approximation  of  I  is  found.  According  to 
[Cl,  pp. 28-38]  the  approximation  of  I,  Iq0  is: 

I00  "  [2’rfA][|Det(Gu)|]'1/2exp(j[kg+(»t/4)(sgn(gulul")  + 

+  s8n<8u2u2"),l)|ui.0>  u2-0  (C2) 


with: 


Det(Cu)  -  gxx  gyy  '  -  Sxy' ’ Bxy' '  “ 

_  *2e  S2B  *2g  S26 

SxSx  SySy  SxSy  SxSy 


(C.  3) 
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avSe 


®ulul  0 

(  cos#  -sin 5 

8xx 

®xy 

cos#  sin# 

0  ®u2u2  *  • 

(  sin#  cos#  4 

■  &xy 

gyy”  . 

-sin#  cos# 

(C . 4a) 

or : 

Sulul"  _  gxx"cos2tf  +  Syy"sln2#  '  2gxy"cosflsinS  (C.4b) 


6u2u2  "  Sxx"sin2#  +  Syy #  ' cos^9  +  2gxy"cos#sin#  (C.4c) 

0  "  (&xx" -6yy"  +  g^'Ccos2#  -  sin2#)  (C.4d) 

Equation  (C.Ud)  yields,  see  also  (Cl,  pp. 28-38]: 

S  -  (l/2)tan'1(2gxy"/(g:)ry"  -  g**")]  (C.5) 

To  apply  the  theory  above,  equation  (52)  is  rewritten: 


E(e) 


£(kx,ky)exp(Jrg(kx,ky))dkxdky 


(C, 6a) 
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with: 


£(kx,ky)  “  £(S) 


(C.6b) 


g(kx,ky)  -  l^x/r  +  kyy/r  +  kzz/r 


(C. 6c) 


The  derivatives  with  respect  to  kx  and  ky  are  found  to  be  (with  kz 
7(k2  -  kx2  -  ky2)  ): 


Skxky"  -  '  [z/r ]  f  (kz2  +  kx2)Az3] 


(C. 7a) 


Ekyky"  -  -[2/r][(kz2  +  ky2)Az3] 


(C.7b) 


Skxky’ '  -  -[z/rj(kxkyAz- 


(C.7c) 


So: 


Det(Gu)  -  Skyky  *  Skxky  ®kxky 

-  [z2/r2]  [k2Az4] 


(C.  8) 


6  -  (l/2)tan’1[2kxky/(ky2  -  l^2)] 


(C.9) 


The  stationary  phase  is  found  for  Vg-O.  This  implicates: 


Vg(kx,ky)  -  fi  ■»  k  -  ISq  -  kftr 


(C . 10a) 
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This  is  true  because: 


?g(kx,ky) 


+ 


(C.lOb) 


or: 


k 


x 


(C . 10c) 


kz  follows  from: 


kz2  -  k2  - 


k  2  -  k  2 
Kx  Ky 


2  2 
,0  o  X*  +  yz 

k  •  **2<— 


(c.iod) 


or: 


kz  -  k7(z2/(x2+y2+z2) )  -  k[z/rj 


(C.lOe) 


Substitution  in  (C.lOe)  yields: 


kx  “  k[x/r] ;  ky  -  k[y/r] 


(C.lOf) 
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so:  ]c  and  £  have  the  same  direction.  This  means  for  (C.8): 

Det(Gu)  -  [z2/r2][k2Az*]  -  (lA^H^A^H^/r2]  - 

-  [lAzO2n(kxO2+ky02+kzO2>/kzO2l[z2/(x2+y2+z2)]  - 

-  [1Az02]i(kxo2+ky02+kz02>/kz02][kz02/(kx02+ky02+kz02)1  " 

-  lAz02  (c-  ID 

Equation  (C.7)  substituted  in  (C.4)  with  x  and  y  replaced  by  k  and  k  , 

x  y 

respectively,  gives: 

Sulul"  “  -[z/(rkz3))[kz2  +  (kxcos9  -  kysin0)2]  (C.12a) 

6u2u2 ' '  "  ' (z/(rk23) ] [kz2  +  (h^cosfi  +  kysin0)2]  (C.12b) 

Because  r  -  J(x2+ y2+z2)  >  0  and  sgn(z)  -  sgn(kz)  -  sgn(kz3),  for  kzeR 
(the  latter  following  from  (C.10)  ): 

sgnlguiux"!  “  s8rU6u2u2">  ”  -1  (C.13) 

and  so  (C.6)  becomes: 

E(E)  ~  [2>r/r]t(K)[  |1Az02|  )  "1/^2exp(  j  [r(k  x/r  +  k  y/r  +  k  z/r)  + 

r  ■>  «  x  y  z. 

+  (x/4) (-1-1)  ]  )  - 

+jkr 

e 

-  - j 2xkcos#t (JO -  (C .  14) 

r 

in  which  use  is  made  of: 


kxOx  +  kyQy  +  kzOz  -  ko'I  "  k*r-rir  -  kr 


(C.13) 
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